DESCRIPTION 



MAGNETIC RESONANCE IMAGING DEVICE, IMAGE DATA CORRECTING 
DEVICE AND IMAGE DATA CORRECTING METHOD 



Technical Field 

[0001] The present invention relates to a magnetic 

resonance imaging device, an image data correcting device 
and an image data correcting method for obtaining a 
magnetic resonance (MR) image of a detected body by 
utilizing a magnetic resonance phenomenon. In particular, 
the present invention relates to a magnetic resonance 
imaging device, an image data correcting device .and an 
image data correcting method for executing a correction of 
collecting data spatially ununiformly deteriorated by a 
nonlinear movement of the detected body when magnetic 
resonance imaging using a spin warp method, a spiral 
method, a radial method, etc. for collecting data filled 
in a k-space by performing a scan of a multi-shot is 
performed. 

Background Art 

[0002] The magnetic resonance imaging is now used as one 
of important medical modalities in a medical spot. A 
system for executing this magnetic resonance imaging is 
called the magnetic resonance imaging device. A high 
frequency magnetic field is applied to the detected body, 
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and the magnetic resonance phenomenon is generated in a 
magnetized spin within the detected body. An echo signal 
generated by this magnetic resonance phenomenon is 
collected. Such operations are executed under a time 
series. An MR image of the detected body is basically 
obtained by processing this echo signal (including 
reconstruction processing) . such a series of operations 
is called a scan. 

[0003] in this magnetic resonance imaging, it is not 

preferable to move the detected body (including a body 
movement due to beating and breathing) during this scan. 
Such a movement becomes a motion artifact and deteriorate 
the quality of the reconstructed MR image. Therefore, the 
restraint of the motion artifact is an important technical 
theme in the magnetic resonance imaging. 

[0004] A method using an inclining magnetic field pulse 
called a navigator is known as one of methods for 
restraining the motion artifact when the magnetic 
resonance imaging is performed by using a pulse sequence 
of a multi-shot type. In this method, the movement of the 
detected body is monitored by using the echo signal (als 
called the navigator echo signal, or briefly called th 
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naviecho) collected by applying this navigator, and the 
collecting data are corrected on the basis of this monitor 
information (see non-patent literatures 1 to 8) . 



[0005] As can be seen from non-patent literatures 1, 2, 
the mode of data for monitoring the movement was first a 
one-dimensional projecting profile. In contrast to this, 
recently, as can be seen from non-patent literatures 4, 5, 
7, a system for two-dimensionally monitoring and 
correcting the movement during the scan is also proposed. 

[0006] However, in the methods reported so far, a rigid 
movement of a head portion, etc., i.e., parallel 
translation and a rotating movement as a rigid body are 
set to objects. With respect to a movement having a 
spatial distribution in the movement, i.e., a nonrigid 
movement, a movement collected within a voxel as in 
diffusion imaging is set to an object (see non-patent 
literature 7) . It is merely reported in application (see 
patent literature 9} able to separate a movement component 
in a time series image of the head portion at a spatial 
frequency. However, it is not reported in a structure in 
which the nonrigid movement for spatially ununiformly 
moving (shifting) each position forming the abdominal part 
as in the movement (motion) of a breathing property of the 
abdominal part, etc. is set to an object. 

[0007] The technical reasons for this unreport will be 

explained. When the rigid body is shifted in a linear 
position in a direction parallel with this inclining 
magnetic field under an inclining magnetic field of 
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constant intensity and constant time in the magnetic 
resonance imaging, this position shift is proportional to 
a shift amount of an average phase in the r ( real ) -space 
(exactly, it is also called a hybrid space (h-space) since 
it is the r-space on only one side axis and is the k-space 
on the other axis) . Accordingly, when the shift is 
nonlinear every position, i.e., the shift amount every 
position is different, it is impossible to know the 
spatial distribution of the movement even when the average 
phase is calculated. 

[0008} On the one hand, the nonlinear movement can be 

corrected if the spatial distribution of the position 
shift caused by the movement can be measured in accordance 
with a changing frequency. It is necessary to newly apply 
a navigator pulse every shot as well as the pulse sequence 
for imaging so as to measure the movement. However, in 
the case of this technique, time-like restriction is added 
to the navi echo collected to calculate a two-dimensional 
distribution of the movement. On the other hand, since it 
is comparatively easy to measure the movement as one 
dimension as projection data of a certain axis direction, 
there are many reports. However, the application of this 
technique is limited to the linear movement as the rigid 
body, and it is difficult to cope with the nonlinear 
movement . 
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[0009] The case for measuring the navi echo as the one- 
dimensional projection data will be described in detail. 
In this case, normally, in addition to the echo for 
imaging, the navi echo of a read-out direction provided 
when a phase encode amount is set to zero every shot and 
the navi echo is collected while applying the inclining 
magnetic field in the read-out direction of the imaging, 
i.e., projection data in the phase encode direction are 
acquired and corrected. In this case, with respect to the 
navi eco of the read-out direction, the projection data 
integrated along the phase encode direction can be 
acquired with respect to each position of the read-out 
direction in timing considered to make the same movement 
as the collecting time of the echo for imaging. Therefore, 
the echo for imaging can be corrected correspondingly to 
each line of each projection data in the read-out 
direction. Accordingly, it is also possible to cope with 
the nonlinear movement in the read-out direction to a 
certain extent (see non-patent literature 1}. 
[0010] However, in the case of the navi echo collection 
of the read-out direction at the imaging time, the 
nonlinear correction with respect to the movement in the 
phase encode direction is impossible. Accordingly, the 
navi echo of the phase encode direction is used. However, 
in this case, differing from the navi echo of the read-out 
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direction, the echo data for imaging acquired in the same 
moving state as the navi echo are data on only a line 
acquired in the same shot as the navi echo. Therefore, in 
the phase encode direction of the echo for imaging, data 
acquired in the moving state different every shot are 
mixed. Therefore, even when these data are corrected to 
data of the r-space (also called the hybrid-space (h- 
space) ) by performing the Fourier transform in the phase 
encode direction as they are, it does not correspond to 
the line of the navi echo Fourier- transformed . 
Accordingly, it is difficult to make the nonlinear 
movement correction of the phase encode direction by the 
navi echo of the phase encode direction. However, in this 
case, the correction can be also made if zero is put close 
to a position except for the shot corresponding to the 
navi echo in the k-space of the echo for imaging, and the 
data transformed to the hybrid space (h-space) are 
corrected every shot and are synthesized. 

[0011] Thus, it is possible to measure the shift of the 
position due to a spatially ununiform movement and the 
distribution of the phase in principle, and make the 
correction in the r-space every voxel. However, a 
measuring technique (a pulse sequence or an external 
monitor) of the movement different from the original 
imaging technique (pulse sequence) is required. 
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Accordingly, it is difficult to make the measurement 
itself in time and technique. In addition to this, it is 
necessary to execute the measurement and the correction by 
a shot number with respect to one image. Therefore, as 
the number of shots is increased, an arithmetic amount 
becomes enormous and nonrealistic. In a case using a one- 
dimensional distribution as projection, the movement of 
the read-out direction of imaging can be also made in 
nonlinearity by measuring the position shift every 
projection line in the r-space and reversely performing 
the shift, but the arithmetic amount is increased, When 
the nonlinear movement of the phase encode direction of 
imaging is set to an object, the arithmetic amount is 
further increased. In this case, no interpolation is 
required in the correction in the k-space. However, in 
the correction in the real space, the interpolation is 
required in the correction of the shift in a degree of 1 
pixel or less . 

[0012] This situation will be described in detail with 

the spin warp method of the multi-shot as an example. In 
the case of the spin warp method of the multi-shot," 
collective data of a certain line number are acquired by 
performing division in the phase encode direction in the 
k-space every one shot. Therefore, data influenced by the 
movement different every shot are mixed in the k-space. 
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In the nonlinear correction of the spatial position, 
turning-back is caused in the image space even when a 

distribution AY(y,n) of the nonlinear shift every position 
y is measured in the r-space every shot, and the corrected 
image is repeatedly made. Therefore, portions of 

different movements are overlapped and cannot be 
distinguished so that no correction can be made in 
principle. With respect to the measurement itself, in the 
case of the nonlinear movement, it is very difficult to 
identify the position relation before and after the 
movement of the same portion every voxel. In particular, 
in the case of the normal spin warp method, data every one 
line are collected every shot and are converted into data 
of the r-space every line, and the shift amount every 
position y is then measured. Thereafter, the position 
correction is made every position y. Accordingly, the 
arithmetic amount becomes enormous in both the measurement 
and the correction. 

[0013] On the other hand, in the k-space, the 

measurement and the correction can be easily made until 
the phase distributions of spatial zeroth and first orders 
The position shift in the r-space becomes a phase shift in 
the k-space. Therefore, it is sufficient to 

arithmetically calculate the product of a phase term at 
the same point of the k-space, and the correction of a 
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subpixel or less in the r-space can be also made. 
Therefore, the arithmetic calculation in the k-space is 
desirable in consideration of convenience of processing. 
However, in the k-space, it is limited to the measurement 
and correction of the linear motion, and the measurement 
and correction of the spatial nonlinear motion are 
difficult . 
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[0014] Disclosure of the Invention 

The present invention is made in consideration of 
the above situation of the prior art, and its object is to 
provide a magnetic resonance imaging device, an image data 
correcting device and an image data correcting method able 
to obtain an MR image by comparatively simply correcting 
echo data spatially ununiformly deteriorated and collected 
by the nonlinear movement of an image pickup part at high 
speed when the image is obtained from data of the k-space 
particularly collected and filled by using the pulse 
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sequence of the multi-shot type in the magnetic resonance 
imaging * 

[0015] To achieve the above object, an image data 

correcting' device in the present invention comprises: 

a movement information acquiring section for 
acquiring movement information showing a spatial 
distribution of the magnitude of a movement in the real 
space of an image pickup part of a detected body; 

a correcting section for making a correction 
different from that of a second area in a first area of 
image data of the image pickup part of the detected body 
collected by a scan of magnetic resonance imaging on the 
basis of the movement information; and 

a synthesizing section for synthesizing respective 
image data of the first area and the second area corrected 
by the correcting section. 

[0016] To achieve the above object, an image data 

correcting device in the present invention also comprises: 

a correcting section for making a correction 
different from that of a second area in a first area of 
image data of an image pickup part of a detected body 
collected by a scan of magnetic resonance imaging on the 
basis of movement information showing a spatial 
distribution of the magnitude of a movement in the real 
space of the image pickup part; and 



a synthesizing section for synthesizing the 
respective image data of the first area and the second 
area corrected by the correcting section. 

[ 0 01 7 ] To achieve the above object, an image data 

correcting device in the present invention also comprises: 

a first data converting section for converting data 
of a first space in a first area and a second area of an 
image pickup part of a detected body into data of a second 
space in a third area and a fourth area; 

a correcting section for making a correction 
different from that of the fourth area with respect to the 
data of the second space in the third area; 

a synthesizing section for synthesizing the data of 
the second space in the third area after the correction, 
and the data of the second space in the fourth area; and 

a second converting section for converting the data 
of the second space after the synthesis into data of the 
first space. 

[0018] To achieve the above object, an image data 

correcting device in the present invention also comprises: 
a first data converting section for converting data 
of at least one area among data of a first space of plural 
areas in an image pickup part of a detected body into data 
of a second space; 



a 



correcting section for correcting the data of the 



second space; and 



a 



second data converting section for converting the 



data of the second space after the correction into data of 
the first space. 

[0019] To achieve the above object, an image data 

correcting method in the present invention comprises: 

a step for acquiring movement information showing a 
spatial -distribution of the magnitude of a movement in the 
real space of an image pickup part of a detected body; 

a step for making a correction different from that 
of a second area in a first area of image data of the 
image pickup part of the detected body collected by a scan 
of magnetic resonance imaging on the basis of the movement 
information; and 

a step for synthesizing the respective corrected 
image data of the first area and the second area. 
[0020] To achieve the above object, an image data 

correcting method in the present invention also comprises: 

a step for making a correction different from that 
of a second area in a first area of image data of an image 
pickup part of a detected body collected by a scan of 
magnetic resonance imaging on the basis of movement 
information showing a spatial distribution of the 
magnitude of a movement in the real space of the imaqe 
pickup part; and 
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a step for synthesizing the respective corrected 



image data of the first area and the second ar 
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t 0021 3 T o achieve the above object, an image data 

correcting method in the present invention also comprises: 

a step for converting data of a first space in a 
first area and a second area of an image pickup part of a 
detected body into data of a second space in a third area 
and a fourth area; 

a step for making a correction different from that 
of the fourth area with respect to the data of the second 
space in the third area; 

a step for synthesi zing . the data of the second space 
in the third area after the correction, and the data of 
the second space in the fourth area; and 

a step for converting the data of the second space 
after the synthesis into data of the first space. 
[0022] To achieve the above object, an image data 

correcting method in the present invention also comprises: 

a step for converting data of at least one area 



among data of a first space of plural areas in an image 
pickup part of a detected body into data of a second 
space; 

a step for correcting the data of the second space; 

and 

a step for converting the data of the second space 



after the correction into data of the first space. 

[0023] To achieve the above object, a magnetic resonance 

imaging device in the present invention comprises: 

a movement information acquiring section for 
acquiring movement information showing a spatial 
distribution of the magnitude of a movement in the real 
space of an image pickup part of a detected body; 

an image data collecting section for collecting 
image data of the image pickup part of the detected body 
by a scan of magnetic resonance imaging; 
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correcting section for making a correction 



different from that of a second area in a first area of 
the collected image data on the basis of the movement 
information; and 

a synthesizing section for synthesizing the 
respective image data of the first area and the second 
area corrected by the correcting section. 

[0024] To achieve the above object, a magnetic resonance 
imaging device in the present invention also comprises: 

an image data collecting section for collecting 
image data of an image pickup part of a detected body by a 
scan of magnetic resonance imaging; 

a correcting section for making a correction 
different from that of a second area in a first area of 
the collected image data on the basis of movement 
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information showing a spatial distribution of the 
magnitude of a movement in the real space of the image 
pickup part; and 

a synthesizing section for synthesizing the 
respective image data of the first area and the second 
area corrected by the correcting section. 

[0025] To achieve the above object, a magnetic resonance 
imaging device in the present invention also comprises:' 

a data collecting section for collecting data of a 
first space in an image pickup part of a detected body by 
a scan of magnetic resonance imaging; 

a first data converting section for converting the 
data of the first space in a first area and a second area 
into data of a second space in a third area and a fourth 
area; 

a correcting section for making a correction 
different from that of the fourth area with respect to the 
data of the second space in the third area; 

a synthesizing section for synthesizing the data of 
the second space in the third area after the correction, 
and the data of the second space in the fourth area; and 

a second converting section for converting the data 
of the second space after the synthesis into data of the 
first space . 

[0026] To achieve the above object, a magnetic resonance 
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imaging device in the present invention also comprises: 

a data collecting section for collecting data of a 

first space in an image pickup part of a detected body by 

a scan of magnetic resonance imaging; 

a first data converting section for converting data 

of at least one area among the data of the first space of 

plural areas into data of a second space; 



a 



correcting section for correcting the data of th 



second space; and 

a second data converting section for converting the 
data of the second space after the correction into data of 
the first space. 

[0027] In accordance with such magnetic resonance 

imaging devices, image data correcting devices and image 



data correcting methods in the present invention, when th 
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image is obtained from data of the k-space particularly 
collected and filled by using the pulse sequence of the 
multi-shot type in the magnetic resonance imaging, echo 
data spatially ununiformly deteriorated and collected by 
the nonlinear movement of the image pickup part at high 
speed are corrected as a sum (synthesis) of data linearly 
corrected in accordance with a moving degree. Thus, an MR 
image comparatively simply corrected at high speed is 
obtained * 

Brief Description of the Drawings 
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[0028] Fig. 1 is a block diagram showing one embodiment 
mode of a magnetic resonance imaging device in the present 
invention and relating to a construction. 

Fig. 2 is a view for explaining a section of the 
abdominal part of a detected body and a movement of each 
part thereof. 

Fig. 3 is a view showing one portion of a pulse 
sequence using a navigator usable in the embodiment mode 
of the present invention* 

Fig. 4 is a view showing one portion of another 
pulse sequence using the navigator usable in the 
embodiment mode of the present invention. 

Fig. 5 is a view for schematically explaining data 
collection, a nonlinear correction and the procedure of a 
reconstruction in accordance with a first embodiment mode. 

Fig. 6 is a graph showing an example of a window 
function W made when data of an area of the k-space are 
divided into three portions. 

Fig. 7 is a view for schematically explaining data 
collection, a nonlinear correction and the procedure of a 
reconstruction in accordance with a second embodiment mode 

Fig. 8 is a view for explaining a coil arrangement 
in accordance with a modified example of the second 
embodiment mode. 

Fig. 9 is a view showing a single coil and a muiti- 
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coil used as a signal receiving RF coil, their sensitivity 
profiles and the magnitude of a non-rigid body motion of 
the detected body in a magnetic resonance imaging device 
in accordance with a third embodiment mode of the present 
invention . 

Fig. 10 is a view showing the window function when 
each - k-space data collected by using two surface coils 
having a sensitivity distribution similar to a one- 
dimensional movement distribution as shown in Fig. 9 are 
divided into two areas so that these data are equivalently 
divided into three portions . 

Fig. 11 is a view for explaining the principle of 
correction processing in each embodiment mode of the 
present invention . 

Fig, 12 is a schematic flow chart showing the 
procedure of processing common to the respective 
embodiment modes of the present invention. 

Fig. 13 is a flow chart showing a procedure when the 
movement of a non-rigid body is three-dimensionally 
corrected by the present invention. 
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ig. 14 is a view showing an example of a moving 



direction seen from the y-direction when there is a three- 
dimensional movement of the non-rigid body in the detected 



body. 



Fig. 15 is a view showing the moving direction seen 



from the x-direct ion of the non-rigid body in the detected 
body shown in Fig. 14. 

Fig. 16 is a view showing a fault image obtained by 
dividing and correcting the k-space data in software using 
the single coil by the present invention when there is a 
movement of the non-rigid body having a one-dimensional 
distribution in a direction perpendicular to a phase 
encode direction* 

Fig. 17 is a view showing a fault image obtained by 
dividing and correcting the k-space data in software using 
the single coil by the present invention when there is a 
movement of the non-rigid body having the one-dimensional 
distribution in a PS direction. 

Fig. 18 is a view showing a fault image obtained by 
dividing and correcting the k-space data by using the 
single coil and including an air portion in software by 
the present invention when there is a movement of the non- 
rigid body having the one-dimensional distribution in the 
PE direction. 

Fig. 19 is a view showing a fault image in which k- 
space data divided into two portions in hardware by using 
the multi-coil having two element coils by the present 
invention are further divided into two portions in 
software by using the window function, and the fault image 
is obtained by linearly correcting the k-space data 
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equivalently divided into three portions when there is a 
movement of the non-rigid body having the one-dimensional 
distribution in the PE direction. 
Best Mode for Carrying Out the Invention 

[0029] The embodiment modes of the present invention 

will next be explained with reference to the drawings. 
[0030] (First embodiment mode) 

One embodiment mode in accordance with a magnetic 
resonance imaging device of the present invention will 
next be explained with reference to Figs. 1 to 5 . 
[0031] Fig. 1 shows the schematic construction of the 

magnetic resonance imaging (MR I) device in this embodiment 
mode * 

[0032] This magnetic resonance imaging device has a bed 

section for arranging a detected body P thereon, a static 
magnetic field generating section for generating a static 
magnetic field, an inclining magnetic field generating 
section for adding position information to the static 
magnetic field, a signal transmi tt ing-rece iving section 
for transmitting and receiving a high frequency signal / 
and a control- arithmetic section for being in charge of 
the control of an entire system and the reconstruction of 
an image , 

[0033] For example, the static magnetic field generating 
section has a magnet 1 of a superconducting system, and a 
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static magnetic field power source 2 for supplying an 
electric current to this magnet 1. A static magnetic 
field H 0 is generated in the axial direction (Z-axis 
direction) of a cylindrical opening portion (a space for 
diagnosis) into which the detected body P is inserted with 
play. An uni llustrated shim coil is arranged in this 
magnet section. In the bed section, a roof plate T 
arranging the detected body P thereon can be inserted into 
the opening portion of the magnet 1 so as to be escaped. 
[0034] The inclining magnetic field generating section 

has an inclining magnetic field coil unit 3 assembled into 
the magnet 1. This inclining magnetic field coil unit 3 
has three sets (kinds) of x, y and z coils 3x to 3z for 
generating the inclining magnetic fields of the X-axis 
direction, the Y-axis direction and the Z-axis direction 
perpendicular to each other. The inclining magnetic field 
section also has an inclining magnetic field power source 
4 for supplying electric currents to the x, y and z coils 
3x to 3z. This inclining magnetic field power source 4 
supplies pulse electric currents for generating the 
inclining magnetic fields to the x, y and z coils 3x to 3z 
under the control of a sequencer 5 described later. 



[0035] The inclining magnetic fields of the direction 
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of the three axes (the X-axis, the Y-axis and the Z-axis) 



as physical axes are synthesized by controlling the puis 
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electric currents supplied from the inclining magnetic 
field power source 4 to the x, y and z coils 3x to 3z. It 
is then possible to arbitrarily set and change logic axis 
directions constructed by a slice direction inclining 
magnetic field G s , a phase encode direction inclining 
magnetic field G E and a read-out direction (frequency 
encode direction) inclining magnetic field G R 
perpendicular to each other. The. respective inclining 
magnetic fields of the slice direction, the phase encode 
direction and the read-out direction are superposed onto 
the static magnetic field H 0 . 

[0036] The signal transmit ting-receiving section has a 

signal transmitting RF (high frequency) coil 7T and a 
signal receiving RF coil 7R arranged in the vicinity of 
the detected body P in a photographing space within the 
magnet 1, and also has a signal transmitter 8T and a 
signal receiver 8R respectively connected to these RF 
coils 7T and 7R. 

[0037] The signal receiving RF coil 7R is a single coil 
as one coil element, or is a multi-coil constructed by 
plural coil elements, and each coil is formed as a surface 
coil. The signal receiving RF coil is arranged along the 
body surface of an abdominal part, etc. as an image pickup 
part of the detected body P (tested person) . For example, 
the signal transmitting RF coil is formed as a coil for 
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the whole body* 

[0038] The signal transmitter 8T and the signal rec eiver 
8R are operated under " the control of the sequencer 5 
described later. The signal transmitter 8T supplies an RF 
(radio frequency) electric current pulse of a Larmor 
frequency for exciting nuclear magnetic resonance (NMR) to 
the signal transmitting RF coil 7T by these operations. 
The signal receiver 8R fetches a magnetic resonance (MR) 
signal {high frequency signal) received by the signal 
receiving RF coil 7R. The signal receiver 8R then 
performs various kinds of signal processings such as pre- 
amplif ication, intermediate frequency conversion, phase 
wave detection, low frequency amplification, filtering, 
etc, with respect to this MR signal. Thereafter, the 
signal receiver 8R performs A/D conversion, and generates 
digital data (original data) of the MR signal. 
[0039] The control-arithmetic section has the sequencer 

(also called a sequence controller) S, a host computer 6, 
an arithmetic unit 10, a memory unit 11, a display device 
12 and an input device 13. The host computer 6 has a 
function for giving a command of pulse sequence 
information to the sequencer 5 by an uniilustrated stored 
software procedure, and generalizing the operation of the 
entire device. 

[0040] As a pulse sequence used in a scan, a pulse 
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series of a two-dimensional or three-dimensional multi- 
shot type using a high speed SE (spin echo) method, a high 
speed EE (field echo) method, an EP1 {echo planar imaging) 
method, a EASE (fast asymmetric (advanced) spin echo) 
method, a GRASE (gradient and spin echo) method, etc* is 
used. Data are divided on the basis of information of the 
r-space even when the correction is executed in the k- 
space. Therefore, the k-space may be sampled in any 
method. Therefore, such a pulse series may be concretely 
a pulse series based on any one of the spin warp method, 
the spiral method and the radial method. Thus, when the 
artifact is localized near a generating source by the 
correction described later, a large correction effect is 
particularly shown. In the case of the magnetic resonance 
imaging in which no movement between generated echoes can 
be neglected, the correction in the present invention can 
be applied even in the scan of a single shot type. 
[0041] The sequencer 5 has a CPU (Central Processing 

Unit). and a memory. The sequencer 5 stores pulse sequence 
information sent from the host computer 6, and controls 
the operations of the inclining magnetic field power 
source 4, the signal transmitter 8T and the signal 
receiver 8R in accordance with this information. Further, 
the sequencer 5 is constructed such that digital data of a 
magnetic resonance signal outputted by the signal receiver 
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8R are once inputted to the sequencer 5, and these digital 
data are transferred to the arithmetic unit 10. Here, the 
pulse sequence information is all information required to 
operate the inclining magnetic field power source 4, the 
signal transmitter 8T and the signal receiver 8R in 
accordance with a series of pulse sequences. For example, 
this pulse sequence information includes information 
relative to the intensities of pulse electric currents 
applied to the x, y and z coils 3x to 3z , application 
times, application timing, etc. 

[0042] The digital data (also called original data or 

raw data) outputted by the signal receiver 8R are inputted 
to the arithmetic unit 10 through the sequencer 5, and are 
arranged in the k-space {also called a Fourier space or a 
frequency space) using its internal memory. These data 
are reconstructed to image data of the real space by 
performing the two-dimensional or three-dimensional 
Fourier transform every one set. 

[0043] The arithmetic unit 10 is also set to be able to 
execute correction processing in the present invention in 
accordance with a predetermined algorithm. This 
correction processing is processing for correcting 
spatially ununiform (nonlinear) deterioration due to a 
movement of the detected body with respect to the image 
data reconstructed at the present time point or already 



acquired as after-treatment. For example, when an image 
of the abdominal part of the detected body P is picked up, 
each sampling position of this abdominal part is spatially 
nonlinearly moved by the movement (body movement) of a 
breathing property. Thus, the acquired image data are 
spatially ununiformly deteriorated. Since this is matters 
forming a central portion of the present invention, these 
matters will be described later in detail* 

[0044] The arithmetic unit 10 can also execute synthesis 
processing and differential arithmetic processing of data 
relative to an image in accordance with necessity- This 
synthesis processing includes processing for making an 
adding calculation every pixel, maximum value projection 
(M1P: Maximum Intensity Projection) processing; etc. 
[0045] The memory unit 11 can store image data performed 
with respect to the above synthesis processing and 
differential processing as well as the reconstructed image 
data. For example, the display device 12 is used to 
display the reconstructed image. Further, parameter 
information desired by an operator, information relative 
to a scan condition, the pulse sequence, the image 
synthesis and the differential arithmetic calculation, etc. 
can be inputted to the host computer 6 through the input 
device 1 3 , 

[0046] Here, the correction processing as after- 
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treatment performed with respect to the above image data 
will be explained. 

[0047] {Nature of nonlinear movement) 

In the spatial movement of a living body as the 
detected body, there are the linear movement and the 
nonlinear movement as already described. "The spatial 



nonlinear movement" means that there is a spatial 
distribution in the amplitude and phase of this movement, 
With respect to a periodic property, there are many cases 
in which the movement due to the breathing property and a 
heartbeat, etc. is periodic. However, a sudden unexpected 
movement (involuntary motion) is non-periodic . 
[0048] When the images of internal organs of the 

abdominal part are formed, each sampling position of an 
image pickup area of the abdominal part is strongly 
influenced by the movement (body movement) of the 
breathing property. This movement of the breathing 
property spatially becomes nonlinear (ununiform) . 
Therefore, in this embodiment mode, when the internal 
organs of .the abdominal part, etc, are diagnosed/ an echo 
signal mixed with a movement component spatially 
ununiformly distributed by such a nonlinear movement of 
the breathing property is collected, A spatially 

ununiformly deteriorated signal component of this echo 
signal is corrected. 
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[0049] The movement can be also classified into rigid 

body deformation causing only parallel translation and 
rotation, and non-rigid body deformation causing linear 
deformation and nonlinear deformation including 
enlargement and reduction and shearing, in accordance 
with the present invention, more accurately, the 
correction can be effectively made with respect to the 
echo signal into which the movement component spatially 
ununiformly distributed by the movement of the non-rigid 
body deformation is mixed. Accordingly, data deteriorated 
by the linear deformation among the non-rigid body 
deformation are set to an object, and the correction can 
be also made. However, in the following description, an 
explanation will be made with the correction about the 
nonlinear deformation as an object. 

[0050] A certain model with respect to the movement of 

the living body is supposed to make this correction (see 
Fig. 2) . In this model, an axial section image of the 
abdominal part is set to an object as a sectional image 
showing the spatially nonlinear movement most preferably, 
but may not be necessarily set. 

[0051] In the axial section image of the abdominal part, 
as in the model shown in Fig. 2, it can be supposed that 
the amplitude of the movement of the breathing property 
can be neglected on the back side, but becomes largest on 
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its opposite abdominal wall side (front side), ' and the 
amplitude therebetween is gradually increased from the 
back side to the abdominal wall side. It is now supposed 
that the amplitude of the movement of a certain two- 
dimensional object body is uniform in the x-direction and 
is increased from the back side to the abdominal wail side 
in the y-direction in proportion to its distance y. When 
it is considered that the axial section shape of the 
abdominal part is approximately elliptical and the 
contrast is approximately symmetrical on the front, rear, 
left-hand and right-hand sides, it is considered that an 
average shift amount of the entire image acquired, in 
certain timing reflects the movement of the central 
portion of a photographed body. At this time, it can be 
considered that the amplitude of the movement in the y- 
direction becomes 2ero at the back side end, and becomes 
maximum and shows a change twice that of the central 
portion at the side end of the abdominal wall. 
[0052] Therefore, the following technique is adopted in 

this embodiment mode. Namely, the k-space is divided by 
using information of the nonlinear movement measured or 
estimated in the r-space, i.e., the magnitude (order) of 
this movement. Phase corrections of a zeroth order 
different from each other are made in these plural divided 
k-spaces . These corrected data are mutually weighted and 
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added* 

[0053] (Detection of nonlinear movement and its 

correction) 

Next, a concrete technique of the detection of the 
above nonlinear movement and its correction will be 
explained. Figs. 3 to 5 schematically show such a 
technique . 

■ [0054] It is now supposed that the image pickup part is 
the abdominal part of the detected body, and a surface 
coil constructed by the single coil is arranged as the RF 
signal receiving coil 7R around this abdominal part (see 
Fig . 5 (a ) ) . 

[0055] As mentioned above, each position of the 

abdominal part is spatially nonlinearly moved in 
accordance with a breath. Accordingly, when the abdominal 
part is diagnosed, it is necessary to spatially perform 
division and processing in accordance with the largeness 
or smallness of the movement. On the other hand, since it 
is desirous to execute the phase correction in the k-space 
from convenience of the processing, it is desirous to 
divide data filled in the k-space dependently on the 
magnitude of the spatial movement. 

[0056] In the data filled in the k-space, position 

information relative to the distribution of the movement 
is already equally scattered. Accordingly, it is 
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difficult to perform the division using the spatial 
information of the k-space itself, 

[0057] Therefore, a method for dividing data of the r- 

space by using the spatial distribution information of the 
movement is adopted. Here, the axial section of the 
abdominal part is set to an object. a) The knowledge 

(presumption) of "the abdominal side is larger than the 
back side' 7 is used as an index of the distribution in the 
monitor of a zeroth order phase and the space dependency 
of the movement. Otherwise/ b) the forward-backward 



direction (direction directed from the back side to th 
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abdominal wall direction) of the detected body is actually 
set to the phase encode direction, and a method for 
monitoring the movement of the abdominal part of this 
direction by a navigator echo {hereinafter called a navi 
echo) responsive to an inclining magnetic field pulse for 
navigation, etc. is adopted. 

[0058] The shift of the average (zeroth order) phase in 
the r-space becomes a phase shift of a first order in the 
k-spa.ce. Further, the shift of the position in the r- 
space becomes a phase shift of the first order in the k- 
space * 

[0059] Accordingly, when a phase error is exp(i(j) 0 ), the 

product of each pixel data and exp(-i<J> 0 ) is arithmetically 
calculated uniformly in both the r-space and the k-space 
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to correct its phase. Further, if the shift distance in 
the r-space is AY in the y-direction, the product with 

respect to exp [ -27cK y AY] is arithmetically calculated in 
the k-space. 

[0060] When a method for monitoring only the average 

amplitude of the movement described in the above item a) 
and linearly externally mounting and correcting an 
amplitude change of the movement in the 'phase encode 
direction is adopted, it is preferable to use a pulse 
sequence to which an inclining magnetic field causing the 
phase shift of a constant value in the phase encode 
direction in the navi echo is applied. Fig. 3 shows an 
example of this pulse sequence. 

[0061] When the projection data of the phase encode 

direction described in the above b) item are collected, a 
similar inclining magnetic field is applied in both the 
phase encode direction and the read-out direction before 
the navi echo is collected. Fig. 4 shows this example. 

[0062] Each of Figs. 3 and 4 shows a pulse series based 
on the high speed SE method to which the spin warp method 
is applied as the pulse sequence. However, only one echo 
with respect to one shot is illustrated from the relation 
of the paper face. 

[0063] The echo signal is collected by using such a 

pulse sequence shown in Fig. 3 or 4, and its echo data are 
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arranged in the two-dimensional k-space {see Fig. 5(b)). 
These echo data are reconstructed to image data of the 
real space performed with respect to the two-dimensional 
Fourier transform (2BFT) by the arithmetic unit 10 (see 
Fig. 5(c) ) . 

[0064] Here, when the data distribution of the k-space 

is divided in accordance with the spatial movement of the 
abdominal part, an example for most simply dividing the k- 
space into two portions in its phase encode direction is 
shown. In this case, the image of the real space is 
multiplied by a window in the y-axis direction 
corresponding to the phase encode direction (more 
particularly, the image is multiplied by a window 
function) . The image is then divided into respective two 
dominant areas in a moving degree in data on the abdominal 
wall side (the front side of the abdominal part) and the 
back side (the rear side of the abdominal part) (see Figs. 
5(c) and 5(d)). The two real spaces divided in these data 
are respectively transformed to two two-dimensional k- 
spaces by two-dimensional inverse Fourier transform 

(2DIFT) {see Figs. 5(e) and 5(f)). 

[0065] On the other hand, the navi echo S navi (K y ,n) is 

collected every shot by executing the above pulse sequence 
(see Fig. 5(g)). Accordingly, the phase error A6 0 {K y ,n) as 
movement information of the abdominal part is measured 
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from this navi echo S navi (K y ,n) (see Fig. 5(h)), 
[0066] This measurement is concretely executed every 

shot n as the following formulas. 
[Formula 1] 

$o (K y/ n) =2tiY (n) K y 

<|>o (K y/ n base ) ^2ttY (n base ) K y 

A(|>o (K yi n} -<j> 0 (K y/ n) -<J> 0 (K y , n base ) 
Here, Y(n) is a shift amount of the movement in the shot n. 
[0067] Next, the phase correction is made with respect 

to only data of the k-space having data of the abdominal 
wall side (see Fig, 5(i)). No phase correction is made 
with respect to data of the k-space having data of the 
back side. Such a phase correction is made by multiplying 
data of each shot every shot n by the following formula. 

[Formula 2] 

exp [-iA<|>o (K y , n) ] 
[0068] The two-dimensional Fourier transform (2DFT) is 

performed with respect to each of these two k-spaces, and 
it is returned to the individual real spatial data. A 
pixel value is mutually added and the final real spatial 
data are obtained (see Figs. 5(j) and 5(k)). 
[0069] Here, the collection of the above echo data 

(including data of the navi echo) and the correction 
processing will be explained by forming numerical formulas. 
In this embodiment mode, this processing is executed by 



the arithmetic unit 10, but one portion of this processing 
may be also entrusted to the host computer 6. 
[0070] Here, the following symbols are set. 
[Formula 3] 

Simg(x / y): echo data of an imaging object, 

Snavi(x,Ky): navi echo of the x-axis (read-out) 

direction performed with respect to one-dimensional 

Fourier transform (FT) , 

Ax navi (x,n): amplitude of movement of the x-axis 

direction sampled in n-th shot 

(x,n) : amplitude of movement of the y-axis 

direction sampled in n-th shot 

<J>Xnavi (Kx, n} : phase of navi echo of the x-axis 
direction sampled in n-th shot 

navi(Ky / n) : phase of navi echo of the y— axis 
direction sampled in n-th shot 

W L ,M r s(y): window function for selecting a control 
area of movement in a human body 

n: shot number (l = l-Nm ax ) 
[0071] First, echo data are collected by using the pulse 
sequence based on the FSE method using the normal spin 
warp method as simplest one echo/shot of the multi-shot 
type shown in the above Fig, 3 or 4, and its image data 
are obtained {processing 1} . Namely, when echo data 
Si,nQ(Kx,Ky) are set, image data P img (x,y) of its real space 
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are obtained as the following formula. 
[Formula 4] 

F i rng ( X , y ) 

= FT2 Dkx^x, Ky-»y [ Simg { Kx , Ky ) ] 
At this time, navi echo S n3Vi (K y ,n) is also simultaneously 
collected* 

[0072] Next, the arithmetic unit 10 selects a spatial 

area in accordance with the order (magnitude) of the 
movement. Namely, a window for dividing the area of the 
k-space is determined (processing 2} * 

[0073] Concretely, both edge positions y m , n and y maK of 

the abdominal part of the detected body in the phase 
encode direction are first detected. If the navi echo is 
collected and the profile of a one-dimensional movement is 

measured and is used, an index AY(y) of the spatial 
distribution of the movement is arithmetically calculated. 
Thus, the profile W{y) of the y-axis direction, i.e., the 
phase encode direction is obtained as the following 
formula . 

[Formula 5] 
W (y) 

=AY(y) / {Ay max -Ay min } 

(or AY may be replaced with phase <j>) 
If necessary, this profile W(y) may be also smoothed and 
set to the profile W(y). 
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[0074] Accordingly, for example, when the k-space data 

are divided into two portions in the phase encode 
direction (here the y-axis direction: see Fig, 2) , the 
window function W LrMrS (y) for performing these two divisions 
becomes the following formulas . 

[Formula 6] 

W L (y)=W(y) 

W s (y)=l-W L (y) 

If the k-space data are divided into three portions {see 
Fig. 6), the window function W L#M , s (y) becomes the following 
formulas . 

[Formula 7] 

W L (y ) =--2 - max [W (y) - 0 . 5, 0] 
W M (y) ={i-w L (y) } (when W(y)>=0.5) 
=W(y) (otherwise) 

w s (y)=l»w L {y)-w M (y) 

[0075] Fig. 6 explains setting of the window function 

^L,n,s(y) when the k-space data are divided into three 
portions. In this figure, as an index showing the 
movement distribution in the y-direction of the r-space, 
its ordinate axis is shown by phase <J) ♦ This ordinate axis 
may be also shown by the position shift AY. 

[0076] Next, when the k-space is divided into three 

portions, the arithmetic unit 10 respectively divides the 
data of the k-space into data 3 L (K x ,K y ), S M (K X/ K V ) and 
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5 S (K X , Ky) of a larger area, a medium area and a smaller 
area of the movement by using the calculated window 
function W L , M , s (y} and. the image data P img (x f y) of the real 
space. Namely, the data of the k-space are divided into 
data of plural areas (processing 3) . In this case, the 
following arithmetic calculations are made. 

[Formula 8] 

Pl (x, y) =W L ( y) P img (x, y) 

(with respect to image data of the larger area of 
movement ) 

Pm (X / y) =W M ( y) P img (x, y ) 

(with respect to image data of the medium area of 
movement ) 

P s (x, y) =W S (y) Pimg (x, y) 

=Pimg (x, y) -P L (x, y) -P M (x, y) 

(with respect to image data of the smaller area of 
movement ) 

S L (K x/ K y } =1 FT x _»kx, y -»Ky [ Pl (x, y ) ] 

Sm (K x/ K y } =1 FT x ->kx, y->Ky [ Pr-i (x 7 y ) ] 

S s (K x/ K y ) =IFT x - >KXf y -*Ky[Ps(x,y) ] 

-S (K x , K Y ) -S L (K x , K y ) -S M (K x/ K y ) 
[0077] Next, the following processings 4 to 6 are 

sequentially repeatedly executed by the arithmetic unit 10 
every shot n, i.e., until shot n=l to N max - 

[0078] First, the phase distribution in the k-space 
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(|>(Ky,n) caused by the movement is calculated from the navi 
echo SnavitKy,!!) {processing 4), 
[Formula 9] 

Pnavi iy, n) = FT K y-*y [ Snavi (K yy 11 } ] 

AY (y, n) -Y ( y, n) -Y {y, n base ) 

AY meari (n) ^JAY (y, n) dy 

AY max (n) =max of [AY(y,n)] 

Here, the echo data of a base are selected on the 
basis of central data (n=N roa .*/2) of the k-space or an 
average of all AY(y,n) . 

[0079] Next, a maximum value (|w(K y/ n) of the phase 

distribution in the k-space is arithmetically calculated 

(processing 5) . Concretely, when a linear model and a 
zero order phase are used, the following formulas are 
formed. 

[Formula 10] 

A<J>o (K y , n) =27cK y AY inean (n) 

A^max (K y/ n) =2 -max of [A<|>o (K y , n> ] 
When the profile AY{y,n) of one-dimensional projection of 
the y-axis direction is used, the following formula is 
formed. 

[Formula 11] 

A<|> max (K y , n) - 2 TcKyAY^ (y, n) 
[0080] Next, the phase correction is made with respect 

to K y in the same shot (processing 6) , Concretely, the 
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following formulas are formed. 
[Formula 12] 

S L .cor ( K x , Ky ) 

=S L (K*, K y ) exp [-iA<j> max (K y/ n) ] 

. cor (K X/ Ky) 

^S H (K x/ K y ) exp [-iA<|> max (K y , n) /2 ] 
Thus, the loop every shot n is terminated. 

[0081] Next, the arithmetic unit 10 synthesizes all the 
divided data of the k- space by the following formula 
(processing 7) . 
[0082] [Formula 13] 

S img .cor ( ^x r ) 

cor ( K x/ K y ) 4-Sm. cor ( K x , K y ) + Ss ( K x , Ky ) 

The correction data S iing , cor <K X , K y ) of the entire k- 
space calculated in this way are reconstructed to image 
data of the real space by the two-dimensional Fourier 
transform (processing 8), If this reconstruction is 
formed by a formula, this formula is provided as follows, 

[Formula 14] 

J- img . cor (x, y) 

— FT2D|<x~>x f Ky-+y t Simg . C or ( K X / Ky ) ] 

[0083] Thus, in accordance with this embodiment mode, 

the k-space data are divided in consideration of a moving 
degree by using information of the movement of the 
abdominal part measured or estimated in the r-space. 
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Further, linear movement corrections different from each 
other according to such a moving degree are executed every 
divided data. The corrected divided data are finally 
synthesized. Here, the linear correction includes all 
corrections defined by affine transformation. 
[0084] Namely, with respect to data deteriorated by the 
nonlinear movement of an image pickup part in which a 
large amount of complicated and arithmetic calculations 
are conventionally required and it is substantially 
difficult to execute these calculations, a nonlinear data 
correction is realized in substitution by dividing the 
image data in accordance with the moving degree, and 
making the linear moving correction with respect to each 
divided data, and synthesizing the divided data after this 
correction. 

[0085] Even when MR imaging of the multi-shot type is 

performed by this substituting nonlinear correction, the 
phase of the data spatially deteriorated by the spatial 
nonlinear movement of the abdominal part, etc. can be 
easily corrected at high speed. As this result, the 
artifact caused by generating the spatial ununiform 
nonlinear position shift and phase shift by an influence 
of such a nonlinear movement can be simply corrected at 
high speed, and image quality can be improved. 
[0086] Here, an explanation of the above pulse sequence 
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of Fig. 3 will be supplemented. This pulse sequence is a 
sequence provided when only a zeroth order (average) phase 
shift is measured in the r-space provided by the movement 
along the phase encode direction of imaging every shot .# 
in the navi echo. An inclining magnetic field (navigator) 
of intensity G m and time T m is applied in the phase encode 
direction Ge before the navi echo S navi (K y/ n} , and an 
average phase (corresponding to a central phase in the k- 
space) in the r-space is measured, and a difference with 
respect to a base shot is calculated. 
[0087] [Formula 15] 

taav±o(n)=arg[ S na vi (K y =0,n) ] 

A(|> n aviO ( Y / n) =<J>naviO (Xl) ~(|>navi0 (n base ) 

[0038] It can be also converted into an average shift 

amount AY 0 (n) in the r-space by using the relation of 
<|>navio (n) = (2TcG m T m ) AYo (n) . Further, A/D sampling with respect 
to the navi echo S na vi(K y/ n) may be narrowly set since it is 
sufficient to be able to measure only the phase of the 
center (DC) of the k-space* 

[0089] In a pulse series portion for the navi echo, the 

FE method able to shorten time may be also used instead of 
the SE method if a signal can be sufficiently secured. 
This is because susceptibility effects caused by the FE 
method and the phase shift originated from the static 
magnetic field are canceled since the difference with 



43 



respect to the echo using the base shot is calculated. 
[0090] (Second embodiment mode) 

Next, one embodiment mode in accordance with the 
magnetic resonance imaging device of the present invention 
will be explained with reference to Figs . 1 and 7. The 
hardware construction of the magnetic resonance imaging 
device in this embodiment mode is the same as the above 
hardware construction described in Fig, 1, and its 
explanation is therefore omitted. 

[0091] In this embodiment mode, a multi-coil is used in 
the signal receiving RF coil 7R in hardware, and an 
efficient correction is realized by linking data 
collection utilizing respective sensitivity areas of two 
coil elements forming this multi-coil, and correction 
processing of the present application invention. When a 
case for picking-up the image of the abdominal part is 
particularly considered and a ghost {a nonlinear component 
caused by the movement) from the abdominal wall side is 
extended onto the back side, the ghost is included in data 
of the lung side. In this case, even when the supposition 
of Fig. 2 with, respect to the movement is correct, the 
fear of whether its influence can be finally neglected is 
left* In this embodiment mode, such a fear can be 
reliably wiped out. 

[0092] As shown in Fig. 7 (a} , in this multi-coil 7R as 



the signal receiving RF coil, two surface coils 7RU, 7RD 
respectively arranged oppositely on the back side and the 
abdominal wall side of the abdominal part are set to 
element coils. A pickup image provided by combining such 
a multi-coil and a dedicated image reconstruction is 
widely practically used in recent magnetic resonance 
imaging to shorten image pickup time and improve SNR 
(signal to noise ratio) . When the image of the abdominal 
part is picked up, there are many cases in which two 
surface coils are arranged on the back side and the 
abdominal wall side as shown in Fig. 7(a). 

[0093] One surface coil 7RU of the abdominal wail side 

has stronger detecting sensitivity on this abdominal wall 
side by arranging the surface coils 7RU , 7RD in this way. 
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herefore, a large weight is applied to the abdominal wall 



side with respect to a signal from the surface coil 7RU, 
and the movement of this abdominal wall side is detected 
as a dominant component. The other surface coil 7RD of 
the back side has stronger detecting sensitivity on this 
back side. Therefore, a component small in the movement 
of the back side is dominantiy detected in a signal from 
this surface coil 7RD . Namely, two sets of k-space data 
divided by reflecting the nonlinear moving degree of the 

respectively collected 
from the beginning by respectively arranging these two 
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surface coils 7RU, 7RL along the abdominal wall side and 
the back side of the abdominal part (see Figs. 7(b), 7(c)) 
[0094] These two sets of k-space data respectively 

correspond to Figs. 5(e), 5(f) in the above first 
embodiment mode. Accordingly, the Fourier transform and 
the window processing corresponding to Figs. 5(c), 5(d) 
are not required. The processings of Figs. 7(d) to 7(h) 
are the same as the above processings of Figs. 7{g) to 
7 (k) , 

[0095] The echo data spatially differently weighted can 

be directly collected from the abdominal part of the 
detected body P by the two surface coils 7RU, 7RD arranged 
in this way. In particular, there are many cases in which 
the surface coils are conventionally arranged on the 
abdominal wall side and the back side in the signal 
receiving RF coil in the diagnosing case of the abdominal 



part, Therefore, the echo data dominantly influenced by 



large movement of the abdominal wall side from the two 
surface coils 7RU, 7RD, and the echo data dominantly 
influenced by a small movement of the back side are 
collected from the beginning, 

[0096] Therefore, the above correction processing is 

performed with respect to only the echo data from the 
abdominal wail side, and the respective k-space data are 
reconstructed and are mutually added. Thus, final image 
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data Pi^g.cor (x, y ) of the real space can be obtained (see' 
Figs . 7 (d) to 7 (h) ) . 

[0097] Accordingly, in accordance with this second 

embodiment mode, differing from the above first embodiment 
mode, even when the ghost is extended by the movement of 
the abdominal wall side until the back side, a generating 
source signal of the ghost is preponderantly included in 
only the echo data from the surface coil 7RU of the 
abdominal wall side. Namely, the influence of the ghost 
included in the echo data from the surface coil 7RD of the 
opposite side is restrained from the beginning. Therefore, 
if the echo data from the surface coil 7RU of the 
abdominal wall side are corrected as they are, it becomes 
equivalent to the correction including the ghost extended 
to the back side far from the abdominal wail side. 
Accordingly, in comparison with the case of the single 
coil adopted in the first embodiment mode, it is possible 
to finally obtain the reconstructing image of image 
quality further improved in processing time and correction 
precision . 

[0098] In the above coil arrangement, as shown in Fig. 8, 
it is also possible to construct an arrangement of four 
element coils in which element coils 7A to 7D are 
respectively divisionally arranged in the vertical 
direction (direction from a head portion to a leg portion) 
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on the back side and the abdominal wall side of the 
abdominal part. In this case, since the movement of the 
upper side (head portion side) of the abdominal part 
corresponding to the vicinity of the diaphragm is large, 
separate corrections are preferably respectively made with 
respect to the four element coils. Thus, the correction 
effect is large with respect to the images of a coronal 
face and- a sagittal face. 

[0099] Further, it is also possible to adopt an 

arranging construction of 1+4 coils in which an 
independent surface coil in accordance with the first 
embodiment mode and the above four element coils are 
combined and arranged. In this case, the movement of the 
upper side (head portion side) of the abdominal part 
corresponding to the vicinity of the diaphragm is large. 
Accordingly, the separate corrections may be also 
respectively made with respect to the four element coils 
of the divisional arrangement. Thus, a large correction 
effect is obtained with respect to the coronal face and 
the sagittal face. 

[0100] Further, the image pickup part may be also set to 
a portion except for the abdominal part. When movement 
components of the image pickup part can be separated to a 
certain extent from the beginning between mutual signals 
outputted from the coil elements forming the multi-coil, 
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the image pickup part except for the abdominal part may be 
also set- For example, when the image pickup part is the 
heart, it is known that the movements in the lower portion 
and the upper portion of the heart are nonlinear by the 
movement of the breathing property even in the same 
heartbeat phase. In this case, plural upper, lower, left 
and right surface coils are arranged so as to surround the 
heart from the body wall to execute the image pickup using 
the multi-coil in accordance with this second embodiment 
mode. Thus,. each surface coil can detect the moving 
degree in a separating state in accordance with its coil 
sensitivity. Accordingly, the technique in accordance 
with this second embodiment mode can be executed. 
[0101] In the pickup image of the abdominal part, the 

effect of the case extending the ghost from the abdominal 
wall side to the back side is also already confirmed by a 
simulation executed by the present inventors. 
[0102] In the measurement of the spatial distribution of 
the movement using the multi-coil, the average shift 
amount of the position may be arithmetically calculated 
from the navi echo outputted by each coil element, and may 
be also used in the correction. In accordance with such a 
construction, the correction can be made even when no 
information relative to the position relation of the coil 
is identified. Thus, generation application can be 
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performed irrespective of parts and coil modes. 
[0103] (Third embodiment mode) 

A magnetic resonance imaging device in accordance 
with a third embodiment mode is executed by combining 
processing based on the multi-coil (i.e., the construction 
explained in the above second embodiment mode) , and the 
division of the moving component using the window 
processing using the image processing (i.e., the 
construction explained in the above first embodiment mode) 
Namely, in the magnetic resonance imaging device in 
accordance with the third embodiment mode, the window 
processing is performed with respect to each image from 
each coil element and data are divided to subdivide the 
moving component ♦ 

[0104] Fig, 9 is a view showing the single ooil and the 

multi-coil used as the signal receiving RF coil, their 
sensitivity profiles and the magnitude of non-rigid body 
motion of the detected body in the magnetic resonance 
imaging device in accordance with the third embodiment 
mode of the present invention. 

[0105] As shown in Fig. 9(a), in the third embodiment 

mode, similar to the second embodiment mode shown in Fig. 
7 fa) , a multi-coil 7RM is used as the signal receiving RF 
coil. As shown in Fig. 9(a), in the multi-coil 7RM, two 
surface coils 7RU, 7RD are set to element coils. The 
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surface coils 77RU, 7RD are respectively oppositely 
arranged on the back side and the abdominal wall side of 
the abdominal part of the detected body P. Further, a 
sleeve-shaped single coil 7RS is also arranged as the 
signal receiving RF coil so as to cover the detected body 
P in accordance with necessity. One or both of the single 
coil 7RS and the multi-coil 7RM can be arbitrarily used as 
the RF coil for signal reception, 

[0106] Fig, 9(b) shows the magnitude of the non-rigid 

body motion in the abdominal part of the detected body P , 
The ordinate axis of Fig. 9(b) shows the position y 

[pixel] of image data shown in Fig. 9(a), and the abscissa 
axis shows the magnitude A(y) [pixel] of the movement in 
the pixel y. 

[0107] As shown by the solid line A(y) within Fig. 9, it 
is supposed that the movement forms a first order 
distribution of the y-direction between y u and y d . Namely, 
it is supposed that the magnitude A(y) of the movement is 
a maximum value Amax in y u , and is 0 in y d . 

[0108] Fig. 9(c) shows respective normalized sensitivity 
distributions Wcu(y), Wcd(y), Wcs (y) of the respective 
surface coils 77RU, 7RD and the single coil 7RS in the y- 
direction. The ordinate axis of Fig. 9(c) shows the 
position y [pixel] of the image data shown in Fig. 9(a), 
and the abscissa axis shows respective normalized 
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sensitivity distributions Wcu(y), Wcd(y), Wcs (y) of the 
respective surface coils 77RU, 7RD and the single coil 7RS 
in the y-direction. 

[0109] As shown in Fig. 9(c), similar to the movement 

distribution, sensitivity distributions Wcu ( y ) , ' Wed ( y ) of 
the respective surface coils 77RU, 7RD are respectively 
set to a first order distribution in the y-direction + The 
sensitivity distribution Wcs(y) of the ■ single coil 7RS is 
constantly set in the y-direct ion . 

[0110] The k-space data collected by such respective 

surface coils 77RU, 7RD or single coil 7RS can be divided 
into plural areas by the window function. In particular, 
the k-space data collected by the two surface coils 77RU, 
7RD are divided into two areas in hardware. Therefore, 
one or both of the respective areas can be further divided 
into plural areas in software by using the window function. 
Further, the k-space data collected by the single coil 7RS 
can be divided into an arbitrary number of areas in 
software by using the window function. For example, the 
k-space data collected by the single coil 7RS can be 
divided into two portions or three portions in software by 
using the window function similar to that of the first 
embodiment mode . 

[0111] Fig. 10 is a view showing the window function 

when each of the k-space data collected by using the two 
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surface coils 77RU, 7RD having sensitivity distributions 
similar to the one-dimensional moving distribution as 
shown in Fig. 9 is divided into two areas so that the k- 
space data are equivalently divided into three portions. 
[0112] In Fig. 10, the axis of ordinate shows the values 
of window functions W L (y), W„(y), W s (y) and the axis of 
abscissa shows pixel positions of the y-direction shown in 
Fig. 9. Two k-space - data collected by the two surface 
coils 77RU, 7RD using the three window functions W L (y), 
W M (y), W s {y) shown in Fig. 10 can be divided into three 
areas constructed by a larger area L, a medium area M and 
a smaller area S of the movement. The correction of 
intensity according to the magnitude of the movement can 
be made with respect to the k-space data of each divided 
area . 

[0113] Namely, the k-space data divided into two 

portions in hardware by the two surface coils 7RU, 7RD are 
further divided into two areas in software by using the 
window function. Further, one portion of the areas can be 
synthesized such that the ' k-space data included in two 
medium areas of the movement among the four areas are 
included in a single area. Thus, one portion of the k- 
space data divided and collected by the plural element 
coils is synthesized with another portion and may also 
form the areas . 
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[0114] It is not limited to the plural k-space data 

divided in hardware by the plural element coils, but one 
portion of the k-space data once divided by the window 
function is similarly synthesized with another portion and 
the correction may be then made. The k-space data can be 
divided into predetermined desirable ununiform plural 
areas by using a more standard window function by 
performing synthesis processing of one portion of the k- 
space data before the correction, and the processing can 
be easily performed, 

[0115] Each window function shown in Fig. 10 can be 

calculated from the moving distribution -W(y) and the 
sensitivity distribution Wc(y) of the surface coils 7RO, 
7RD respectively normalized by the following formulas. 
[0116] [Formula 16] 
W L (y) =Wc (y) *W (y) 

= [W (y) ] 2 (if Wc (y)=W (y) } 

W M (y)=Wc (y) * [1-W (y) ]+ [l-Wc (y) ] *W (y) 

[W (y) * { 1-W (y) } ] (if Wc(y)sW(y)) 

w s (y) -i~w L (y) -w M { y) 

-{1-W{y} } 2 (if Wc{y)=W(y) ) 

[0117] Similar to the case for dividing the k-space data 
into three portions in software in the first embodiment 
mode, respective image data Pimg.u(x,y), Pimg*d(x,y) of 
the real space obtained by the two surface coils 77RU, 7RD 
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are multiplied by the above window functions W L (y), W M {y), 
W s (y), and three image data P L (x,y), P M (x,y), P s (x, y) of 
the larger area, the medium area and the smaller area of 
the movement are generated. The two-dimensional inverse 
Fourier transform is performed with respect to each 
generated image data so that k-space data 3 L (Kx,Ky)., 
S H (Kx,Ky), S S (KX/Ky) divided into three areas in accordance 
with the magnitude of the movement can be obtained. 
[0118] Further, with respect to the respective k~space 

data S L (Kx,Ky), S^(Kx,Ky), S s (Kx,Ky), the phase correction 
is executed at different intensities in accordance with 
the magnitude of the movement. The two-dimensional 

Fourier transform is performed with respect to k-space 
data Sa. mg . cor (Kx, Ky) provided by synthesizing the respective 
k-space data S L . cor (Kx, Ky ) , S H . C or (Kx, Ky } , S s (Kx,Ky) after the 
phase correction. Thus, image data P ±mg . cor (x, y } of the real 
space after the correction can be obtained. 

[0119] Thus, in accordance with the third embodiment 

mode, the k-space data can be divided in hardware and 
software by using the plural element coils and the window 
function. Accordingly, both merits obtained by the first 
embodiment mode and the second embodiment mode can be 
given. Further, the k-space data are easily divided into 
predetermined desirable areas by simple processing and the 
correction of intensity according to the magnitude of the 
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movement can be then made by synthesizing one portion of 
the k-space data once divided in hardware and software. 
[0120] A common principle of the correction processing 

of data in each embodiment mode mentioned above will be 
generalized and explained. Fig. 11 shows the principle of 
the correction processing explained in each of the above 
embodiment modes. 

[0121] First, information Input 1 of the movement of the 
non-rigid body is acquired by a certain method. For 
example, the movement may be measured from the navi echo, 
and the information Input 1 of the movement may be also 
obtained by forming a model. Next, image data Input 2 
spatially deteriorated by the movement are spatially 
divided into plural components in accordance with the 
magnitude of the movement by using the information Input 1 
of the movement. Thus, plural k-space data respectively 
divided from each image data are obtained. This division 

(processing 1) of the k-space data may be performed in 
software by the window function, and may be also performed 
in hardware by using the multi-coil. Further, the k-soac^ 
data divided in hardware by using the multi-coil may be 
also further divided in software by the window function. 

[0122] Next, with respect to each component #i(l<i<N) of 
the k-space data divided into N-portions, a linear 
correction is made by linear transformation according to 



56 



the magnitude of the movement (processing 2} . Affine 
transformation (parallel translation) , a rotating movement 
(rotation) , enlargement and reduction (scaling) and 
shearing deformation (shear) are included in this linear 



ransformation. Namely, the linear correction of stronger 



intensity is made with respect to a larger component of 
the movement. In contrast to this, the linear correction 
of weaker intensity is made with respect to a smaller 
component of the movement. With respect to the small 
component of the movement, a correction for setting the 
intensity to zero, i.e., a non-correction may be also set. 
Next, synthesis processing (processing 3) of each 
component of the k-space data after the correction is 
performed. Thereafter, image data Output after the entire 
correction are obtained by FT processing from the k-spac.e 
data after the. synthesis . 

[0123] The flow of the correction processing explained 

in each of the above embodiment modes is generalized and 
is shown in Fig, 12. 

[0124] First, an image pickup part is scanned and imaged 
by the magnetic resonance imaging device/ and echo data 
are collected (step SI) . Navi echo data are collected 
together with this collection of the echo data, or are 
independently collected (step SI')* Data already 

collected may be also adopted in these echo data and navi 



echo dat 

[0125] Next, such processing is entrusted to the 

arithmetic unit 10 of the magnetic resonance imaging 
device. Namely, an image reconstruction for also 

performing the correction processing according to the 
nonlinear movement of the abdominal part, etc. in the 
present invention is executed by the arithmetic unit 10 on 
the basis of various kinds of control information (a 
dividing area number of the movement, the existence of 
window processing, the image pickup part, a body position, 



the kind of coil (the single coil or the multi-coil) ; 



a 



coil arranging position, the- kind of pulse sequence, etc.) 
[0126] First, the arithmetic unit 10 judges whether 

software data division according to the moving degree is 



required or not (step S2). This judgment is made i 



n 



consideration of the kind of a used signal receiving RF 
coil, the coil arranging position, etc, 

[0127] When this judgment is YES, the two-dimensional 

Fourier transform {step S3; see processing c of Fig. 5), 
the window processing (step 54: see processing d of Fig. 
5) , the two-dimensional inverse Fourier transform 
processing (step S5: see processings e, f of Fig. 5), the 
detection of movement information in the r-space (step S6: 
see processings g, h of Fig. 5), the correction of 



deterioration due to the nonlinear movement of data of t! 



ne 



b8 



k-space (step 57: see processing i of Fig. 5), and 'image 
generation using the two-dimensional Fourier and addition 
processing (step S8 : see processings j, k of Fig, 5) 
explained in the above first embodiment mode are executed. 
[0128] In contrast to this, when the judgment in the 

step S2 is NO, processing using the multi-coil explained 
in the above second embodiment mode as the signal 



receiving RF coil is executed. Namely, the processings of 
steps S3 to S5 are omitted, and the detection of movement 
information in the r-space (step S6: see processings d, e 
of Fig. 7) , the correction of deterioration due to the 
nonlinear movement of data of the k-space (step S7; see 
processing f of Fig. 7), and image generation using the 
two-dimensional Fourier and addition processing (step 38 : 
see processings g, h of Fig, 7) are directly executed. 
[0129] Collection data spatially influenced by the 

nonlinear movement, or data spatially directly reflecting 
the nonlinear movement are collected along such a flow. A 
nonlinear component is replaced with a sum of linear 



omponents. It is then possible to reconstruct an MR 



image from which the influence of the nonlinear movement 
of the image pickup part is reliably removed at high speed 
by a small number of comparatively simple processings. 
[0130] In recent years, as can be seen in the above non- 
patent literature 8 y the following technique is reported 
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in the field of the magnetic resonance imaging with 
respect to the heart. Namely, in this technique, a three- 
dimensional position of the diaphragm is calculated in 
real time from the navi echo monitored just before, and 
inclining magnetic field intensity and an RF frequency are 
controlled so that a slab face three-dimensionally follows 
in conformity with the movement of the heart, and the 
motion artifact and defocusing are reduced. This 
technique is an excellent method since no after-treatment 

■ 

is required and there is no image deterioration if it is 
limited to the affine transformation (parallel 
translation) as linear transformation, the enlargement and 
reduction (scaling), and the enlargement and reduction of 
the shearing deformation (shear) , 

[0131] When this technique is compared with the 

technique of the present invention, the shift due to the 
nonlinear movement can be first corrected as well as the 
linear transformation in the technique of the present 
invention. Namely, in the position shift due to the 
movement, there is the predominance of being also 
applicable to the spatial distribution of a high order 
unable in the linear transformation. Further, in the 
technique of the present invention, the division is 
performed after the image reconstruction before the 
correction even when the window processing is required. 
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Accordingly , the arithmetic calculation becomes simple 
although it is after-treatment. Further, after the 

collection is terminated, a corrected image can be 
instantly outputted. Further, when the multi-coil is used, 
the arithmetic calculation is further simplified. If the 
navigation of the movement in a shot unit is executed just 
before the present imaging, there is also a predominance 
property able to perform the correction processing in real 
time . 

[0132] In each of the above embodiment modes, the 

explanation has been made with respect to the one- 
dimensional correction processing, i.e., a case in which 
the moving direction is only the y-direction and amplitude 
shows a one-dimensional distribution. However, the 

measurement (detection) of movement information and its 
correction processing may be also thr ee-dimensionally 
performed in each of the x, y and z directions. In 
particular, the movement information can be measured on 
the basis of the phase distribution in the k-space or the 
shift of the position of at least one direction in the r- 
space . 

[0133] Here, a case for three-dimensionally correcting 

the movement of the non-rigid body will be explained, 
[0134] When the motion of the non-rigid body is 

generalized and both the direction and the magnitude of 
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the shift due to the movement are different every place 
and every time, i.e., when the direction and displacement 
of the shift are a four-dimensional distribution, it is 
sufficient to set a weight function for dividing data to a 
high order in accordance with the distribution of the 
shift due to the movement- Concretely, when the 

displacing direction is mixed in the three directions AX, 
AY, AZ, the window functions are separately calculated 
with respect to the three directions in accordance with 
the magnitudes of the separate independent movements of 
the three directions constructed by the x, y and z 
directions measured or presumed every shot. Thus, the 
division is performed in accordance with the magnitude of 
the movement, and the correction is additionally executed. 
[0135] For example, the window function is made every 

magnitude of the movement (L, M, S in the three divisions) 
every direction of the movement (x, y, z in the 3D space) . 
After the correction is made every component divided by 
the magnitude of the movement with respect to each 
direction, the correction of entire data can be made by 
making an additional calculation. An algorithm example of 
this case is shown below. 

[0136] Fig. 13 is a flow chart showing the procedure of 
a case for three-dimensional ly correcting the movement of 
the non-rigid body by the present invention. In Fig. 13, 
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a symbol provided by adding a numeral to S shows each step 
of the flow chart. Fig, 14 is a view showing an example 
of a moving direction seen from the y-direction when there 
is a three-dimensional movement of the non-rigid body in 
the detected body P. Fig. 15 is a view showing the moving 
direction seen from the x-direction of the non-rigid body 
in the detected body P shown in Fig. 14. 

[0137] First, in a step S10, ' spatial distributions 

AX{x,y,z,n), AY{x,y, z,n), AZ(x,y,z,n) every three 
directions of a displacement vector (X,Y f Z) (see Fig. 14} 
from a certain reference every shot n in the shift of a 
certain position (x,y,z) due to the movement are obtained. 
[0138] Next, in a step 311, maximum displacements 

AX max (x, y, z) , AY ma * (x, y, z ) , AZ max {x, y , z ) in the three 
directions in the 3D space are calculated by the following 
formulas . 

[0139] [Formula 17] 

AX max (x, y, z) =max [ Ax (x, y, z r n) 3 

AY max (x, y, z) -max [AY (x, y, z, n) ] 

AZ max (x, y, z ) =max [ Az (x, y, z, n} ] 
[0140] Next/ in a step S12, respective 3D division, 

window functions (weight functions) Wx(x,y,z), Wy(x,y, z), 
Wz(x,y,z) of the three directions are calculated by the 
f ollowina formulas . 
[0141] [Formula 18] 
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Wx (x,y, z) =AX (x, y, z) /AX max (x, y, z) 

Wy (x, y, z ) ==AY (x, y, z) /AY max (x, y, z) 

Wz (x, y, z } =AZ (x, y, z) /AZ inax (x, y, z) 
[0142] Next, in a step S13, image data P(x,y,z) of the 
real space are divided by the window function of each 
direction by the following formulas with respect to the 
three direction components. Indices of L (large), M 
(medium) and S (small) are set to be added every magnitude 
of the movement. 
[0143] [Formula 19] 

Px L ,M, 2 (x,y,z)=P(x,y, z)Wx L , M(S (x, y, z) 

Py L ,M,s (x, y) =P (x, y, z) Wy L , M ,s (x, y, z) 

Pz l ,m,s (x, y) =P (x, y, z) Wz L , M , s (x, y, z) 
[0144] Next, in a step S14, as shown in the following 

formulas, the inverse Fourier transform is performed with 
respect to the divided image data Px L , M , s (x, y, z ) , 
Pvl,m,s (x, y) , Pz L , H , s (x,y) so that k-space data Sx L , M;S (x, y, z) , 
SyL,M,s (x, y) , Sz L , M ,s(x,y) are obtained. 
[0145] [Formula 20] 

Sx L ,M,s (k x , k y , k z ) =IFT [Px l ,m, s (x, y, z) ] 

Sy L ,M,s (k x , k y , k z }=IFT[Py L , M/S (x,y, z) ] 

Sx L ,„, s (k x , k y , k z )-IFT[Py L , MfS (x,y, z) ] 
[0146] Next, in a step S15, the phase correction is made 
with respect to the k-space data every three components as 
shown in the following formulas. 
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[0147] [Formula 21] 
Correction of movement maximum portion 
Sx.cor . L , (k x , k y , k z } =Sx. L (k x , k y/ k z ) -exp [-ic|)x maJ: (k x , k y , k z ) ] 
Sy.cor. L , (k x , k y , k z ) =Sy . L (k x , k y , k z } * exp [ - i(|>y max ( k x , k y ,k z ) ] 
Sz .cor . L (k x , k y/ k z ) =Sz. L (k x , k y , k z ) ,exp [ -i<j>z max ( k x , k y , k z ) ] 
Correction of movement medium portion 

Sx .cor .„, (k x , k y/ k z ) =Sx. M (k x , k y/ k 2 ) * exp [ -i4tx max (k x , k y , k z ) /2] 
Sy.cor.„, (k x , k y , k z ) =Sy .„ ( k x , k y/ k z ) *exp [~i<|>y max { k x/ k y , k z ) /2] 
Sz.cor. H (k x , k y ,k z )=Sz. M (k K , k y , k z ) *exp [-icj>z max ( k x , k y , k z ) / 2] 
In this case, the following formulas are formed. 
<f> x max (k x/ k y , k z ) =27rk x AX max (x, y, z) /K x 
4>ymax (k x/ k y , k z ) =27rk y AY Inax (x, y, z) /K y 
<J>Zma X (k x , k y , k z ) =27tk z AZ max (x, y, z ) /K z 
[0148] Next, in a step S16, as shown in the following 

formulas, correction components Sx . cor . L , ( k x , k y , k z ) , 
Sy.cor. L , (k x , k y , k z ) , Sz .cor . L (k x , k y , k z ) , Sx.cor. M , (k x , k y , k z ) , 
Sy . cor . m, ( k x , k y , k z ) , Sz . cor . M ( k x , k y , k z ) and non-correction 
components Sx . s ( k x , k y , k z ) , Sy . s ( k x , k y , k z } , s z . s ( k x , k y , k z ) of 
the k-space data are synthesized. Concretely, the 

respective movement components in the correction 
components are added every three directions, and 
components Sx.cor, Sy.cor, Sz.cor of the three directions 
after the addition are then further added. 
[01491 [Formula 22] 
Sx .cor=Sx. cor . L , {k x , k y , k z ) fSx.cor.„, (k x/ k y , k z ) +Sx. s (k x , k y , k z ) 
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Sy.cor-Sy.cor . L , (k x , k y , k 2 ) +Sy.cor ." M| ( k x/ k y , k z ) +Sy . s ( k x , k v/ k z ) 
Sz. cor-Sz. cor . Lf {k x , k y , k z ) +Sz. cor . M , (k x , k y , k z ) +Sz . s {k x/ k y , k x ) 
S „ cor-Sx . cor+Sy . cor+Sz . cor 

Thus, the k-space data S.cor three-dimensionally 
corrected can be obtained. 

[0150] Further, in a step 517, it is possible to obtain 

3D image data P. cor in which the correction of the three- 
dimensional movement is made by 3DFT of the k-space data 
after the synthesis* 

[0151] In each of the above embodiment modes,, the 

technique for detecting the movement information is not 
limited to a technique using the navi echo collected by 
applying the inclining magnetic field as a navigation 
pulse, but may be also set to a technique using external 
optical, magnetic and mechanical sensors. In this case, 
even when no measurement can be made by such a sensor 
until the spatial distribution of the movement, the 
movement of only the abdominal wall is monitored, and the 



spatial distribution may be also presumed in combinati 



on 



with a modeling result of a body movement distribution 
measured or presumed in advance. Thus, since the 

measurement parallel to the data collection of the 
magnetic resonance imaging can be made, it is also 
possible to cope with high speed imaging of the heart, etc. 



[0152] Further, it is also possible to perform th 



e 
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execution by a combination in which control on the fly of 
high frequency, the inclining magnetic field, etc. is 
performed with respect to only a component able to be 
corrected by the linear transformation, and only the 
remaining nonlinear components of the second order or more 
are corrected by the technique in the present invention. 
[0153] Further, in the above correction processing, when 
a measured phase error exceeds a predetermined limit value 
able to be coped as the device, threshold value processing 
may be also performed at this limit value* 

[0154] On the other hand, in each of the above 

embodiment modes, when the image data to be corrected, or 
required image data are one portion of a divided area, the 
divided image data after the correction are not 
necessarily synthesized. 

[0155] In addition, in each of the above embodiment 

modes, after the image data in the r-space are divided, 
plural k-space data are generated by converting each 
divided image data. However, data of a separate space may 



be generated from the data of the r-space by conversi 



on 



able to take over the nature of the data of the r-space, 
and the k-space data may be also generated by converting 
the data of this separate space. Further, conversely, the 
conversion is once performed from the data of the r-space 
to the data of another space obtained by conversion able 
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to take over the nature of the data of the k-space, and 
the k-space data may be also generated from the data of 
this another space. 

[0156] Typical effects common to the above embodiment 

modes and modified examples can be summarized as follows. 
[0157] If the spatial distribution of the one- 

dimensional (projection) or two-dimensional/ three- 

dimensional movement (the position shift or the phase 
shift) can be measured, the phase or a signal value of 
data deteriorated by the nonlinear movement can be 
comparatively simply corrected by arithmetically 
calculating a linear sum of at least two different 
correction data. In the case of the correction of the 
position shift, processing can be performed at high speed 
since the measurement and the phase correction in the same 
wave number in the k-space are made in comparison with a 
case in which the measurement of the position shift and 
the correction processing are performed by crossing a 
voxel in the r-space. The correction in the r-space is 
suitable in the case of the phase correction within the 
voxel. However, in comparison with the case processed 
every voxel, the arithmetic calculation is simple and can 
be made at high speed, 

[0158] Further, even when no spatial distribution of the 
movement is strictly measured, the entire spatial 
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distribution is presumed as a model even in the nonlinear 
case from a whole average and measurement data of the 
movement of a typical portion such as the abdominal wall, 
etc. Accordingly, when the measurement is made by the 
navi echo for simplifying a measuring system, the sequence 
is shortened and a high speed operation can be performed* 
For example, if it is an axial section of the abdominal 
part, it can be supposed that the position shift due to 
the movement is gradually increased from the back to the 
abdominal wall side. 

[0159] Further, even when the single coil is used, it is 
sufficient to reconstruct data by a division number of the 
k-space data. Accordingly, the high speed operation can 
be performed in comparison with the correction in the r~ 
space. In addition, when the image of the abdominal part 
is picked up by using the multi-coil, the echo data can be 
directly acquired. Accordingly, correction accuracy is 
improved, and the high speed operation can be further 
performed since no window division is required* 
[0160] Here, a fault image in the abdominal part of the 
detected body P obtained by making the correction by the 
present invention is shown in comparison with a 
conventional fault image obtained without making the 
orrection . 

[0161] Fig. 16 is a view showing a fault image obtained 



by dividing and correcting the k~space data in software by 
using the single coil by the present invention when there 
is a movement of the non-rigid body having a one- 
dimensional distribution in a direction perpendicular to 
the phase encode (PE) direction. 

[0162] Fig. 16(a) shows an image obtained without making 
the correction. Fig. 16(b) shows an image obtained by a 
conventional uniform mean shift correction. Fig. 16(c) 
shows an image obtained by respectively linearly 
correcting each k-space data divided into two portions in 
software by the present invention. Fig. 16(d) shows an 
image obtained by respectively linearly correcting each k- 
space data divided into three portions in software by the 
present invention * 

[0163] When the conventional uniform mean shift 

correction is made with respect to the image before the 
correction shown in Fig. 16(a), it can be the that the 
correction is insufficient although the ghost of the 
central portion in the forward-backward direction of the 
detected body P conformed to a correcting amount is 
reduced as shown in Fig. 16(b) * 

[0164] On the other hand, as shown in Fig. 16(c), when 

data are divided into two portions in software by the 
present invention and are respectively linearly corrected 
at different intensities, a reduction of SNR and 
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correction insufficiency of the ghost are seen in an 
intermediate (intermediate in the forward-backward 
direction of the detected body P) portion of a window. 
However, it can be confirmed that the ghost is notably 
reduced in comparison with the image obtained by the 
conventional uniform mean shift correction* Further, as 
shown in Fig. 16{d), when data are divided into three 
portions in software by the present invention and are 
respectively linearly corrected at different intensities, 
the reduction of SNR and the ghost in the intermediate 
portion in the forward-backward direction of the detected 
body P are sufficiently corrected, and a n-otable 
correction effect can be confirmed, 

[0165] Fig. 1.7 is a view showing a fault image obtained 

by dividing and correcting the k-space data in software by 
using the single coil by the present invention when there 
is a movement of the non-rigid body having a one- 
dimensional distribution in the PE direction. 
[0166] Fig. 17(a) shows an image obtained without making 
the correction. Fig, 17(b) shows an image obtained by the 
conventional uniform mean shift correction. Fig. 17(c) 
shows an image obtained by respectively linearly 
correcting each k-space data divided into two portions in 
software by the present invention. Fig, 17(d) shows an 
image obtained by respectively linearly correcting each k- 
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space data divided into three portions in software by the 
present invention. Fig, 17(e) is a view showing- a one- 
dimensional distribution of the movement {ordinate axis: 
pixel position, and abscissa axis: normalized magnitude of 
movement ) . 

[0161] When the conventional uniform mean shift 

correction is made with respect to the image before the 
correction shown in Fig. 17(a), it- is seen that the ghost 
including a newly generated ghost is generated before and 
after in the same degree as shown in Fig. 17(b). 
[0168] On the other hand, as shown in Fig. 17(c), it can 
be confirmed that the ghost is reduced when data are 
divided into two portions in software by . the present 
invention, and are linearly corrected at different 
intensities. Further, as shown in Fig. 17(d), when data 
are divided into three portions in software by the present 
invention and are linearly corrected at different 
intensities, the ghost is further reduced. It can be 
further confirmed that the reduction of a signal in an 
intermediate portion of the abdominal wall and the back is 
improved. 

[0169] Fig. 18 is a view showing a fault image obtained 
by dividing and correcting the k-space data by including 
an air portion in software using the single coil by the 
present invention when there is a movement of the non- 
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rigid body having a one-dimens ionai distribution in the PE 
direction , 

[0110] Fig. 18(a) shows an image obtained by 

respectively linearly correcting each k-space data divided 
into two portions by including the air portion in software 
by the present invention. Fig. 18(b) shows an image 
obtained by respectively linearly correcting each k-space 
data divided into three portions by including the air 
portion in software by the present invention. 
[0171] In the air portion, the ghost flying from a 

moving portion is dominant. Therefore, when the image 
data are corrected ■ as maximum (L) data in the movement by 
including the air portion, the ghost is reduced as shown 
in Figs. 18(a) and 18(b) in comparison with Figs. 17(c) 
and 17 (d) including no air portion, 

[0172] Fig. 19 is a view showing a fault image when 

there is a movement of the non-rigid body having a one- 
dimensional distribution in the PE direction. The k-space 
data divided into two portions in hardware by using the 
multi-coil having two element coils by the present 
invention are further divided into two portions in 
software by using the window function. The fault image is 
obtained by linearly correcting the k-space data 
equivalently divided into three portions* 

[0173] Fig. 19(a) shows an image obtained by making the 
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linear correction without dividing the k-space dat 



divided into two portions in hardware by the two element 
coils in software. Fig. 19(b) shows an image in which the 
k-space data divided into two portions in hardware by the 
two element coils are further divided into two portions in 
software, and the image is obtained by correcting the k- 
space data equivalentiy divided into three portions. 
[0174] Fig. 19{a) shows the image obtained by correcting 
only a dominant data side of the movement among the data 
obtained by the two element coils. However, in comparison 
with the image of Fig, 17(c) obtained by correcting data 
divided into two portions in software by using the single 
coil, the ghost of the back side spatially distantly 
flying is reduced. Further, Fig. 19{b) shows the image 
obtained by correcting data equivalentiy divided into 
three portions by dividing data obtained by the two 
element coils into two portions in software. In 
accordance with Fig. 19(b), the reduction of SNR of the 
intermediate portion is further improved. It can be also 
confirmed that a most preferable correcting effect is 
obtained in comparison with other images shown in Figs. 
16(a), (b), (c), (d), Figs. 17(a), (b) , (c), (d) , Figs, 
18 (a) , (b) and Fig. 19(a) . 

[0175] In corrections with respect to an image before 

the correction, an image obtained by the uniform mean 
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shif t correction/ and . an image in a case of the movement 
one-dimensionally distributed in the read-out (RO) 
direction perpendicular to the PE direction, there is no 
difference between the case using the single coil and the 
case using the multi-coil, 

[0176] The present invention is not limited to the 

descriptions of the above embodiment modes, but can be 
embodied by modifying constructional elements in the scope 
not departing from its gist at an embodiment stage. 
Further, various kinds of inventions can be formed by a 
suitable combination of plural constructional elements 
disclosed in the above embodiment modes. For example, 
some constructional elements may be deleted from all the 
constructional elements shown in the embodiment modes* 
Further, constructional elements over different embodiment 
modes may be also suitably combined. When the image of 
the heart is formed, this image is influenced by both the 
movement provided by a heartbeat and the movement of the 
breathing property. However, a method for monitoring a 
breathing period by the navi echo and dynamically 
controlling the position of a slice face in addition to 
heartbeat gating is adopted. 



75 



